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Abstract. We propose a bridge between the theory of exactly solvable models 
and the investigation of traffic flow. By choosing the activities in an apropriate 
way the dimer configurations of the Kasteleyn model on a hexagonal lattice can 
be interpreted as space-time trajectories of cars. This then allows for a calculation 
of the fiow-density relationship (fundamental diagram). We further introduce a 
closely-related cellular automaton model. This model can be viewed as a variant of 
the Nagel-Schreckenberg model in which the cars do not have a velocity memory. It 
is also exactly solvable and the fundamental diagram is calculated. 
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1 Introduction 



In the past years the investigation of traffic flow using cellular automata, has become 
quite popular (see e.g. and references therein). In contrast to the hydrodynam- 
ical description p|, Q this approach is microscopic in the sense that individual 
cars can be distinguished similar to the Follow-the-Leader models [|, 0]. The 
simplicity of this approach allows for a number of new applications which cannot 
be incorporated easily in the more complex theories of traffic dynamics. Especially, 
large scale computer simulations now are possible [§, Q which are effective even in 
the case of complicated street networks. 

In this paper we propose a different approach which uses methods from statistical 
mechanics. It is well known that a number of exactly solvable models of statistical 
mechanics have a graphical interpretation in terms of (closed) graphs. We suggest 
to interpret such graphs as trajectories of cars in traffic. Using the known partition 
function of an exactly solvable free-fermion model, the Kasteleyn dimer model on a 
hexagonal lattice, we then are able to calculate the properties of the corresponding 
traffic model. In this way we not only obtain informations about the fundamental 
diagram (density vs. flow) but also about the correlation functions. Finally, we 
discuss the relation of this model to the approach using cellular automata. 



2 Kasteleyn model 

We suggest and study here a simple statistical two-dimensional lattice model, the 
configurations of which can be mapped onto the trajectories of a discrete (stochas- 
tic) one-dimensional traffic problem on a ring ("Indianapolis situation"), see Fig. 
1. Our traffic model has the following properties: (i) the cars do not collide, the 
minimum distance between any two cars being not less than one lattice spacing; 
(ii) the velocities of the cars are (in principle) unlimited. The principle difference 
between the cellular automaton models and the present approach is that the sta- 
tistical weights are ascribed globally to each allowed set of trajectories, rather than 
determined locally by the transition probabilities from one space configuration to 
another. 

The statistical model is actually an interpretation of the dimer model on the 
hexagonal lattice with anisotropic activities, which has been suggested by Kasteleyn 



|Q] and studied in a series of papers, see 11| and references therein. The hexagonal 



lattice we represent as a decorated square lattice with two sites in each vertex shown 
by a circle in figure la. We consider the horizontal direction of the square lattice as 
the space axis and the vertical one downwards as the time axis of the traffic problem 
as usual in traffic theory. The above interpretation allows us to define an one-to-one 
mapping of each dimer configuration on the decorated lattice. Fig. la, onto a set of 
trajectories on the square lattice. Fig. lb. Under that mapping the circles of the 
decorated lattice are mapped onto the sites of the square lattice; a horizontal dimer 
on the decorated lattice corresponds to a unit space step of some trajectory, and a 
vertical dimer corresponds to a unit time step. The fact that the Kasteleyn model 
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belongs to the class of free-fermion models ensures continuity and non-intersection 
of the trajectories. The number of trajectories and the average velocity of the cars 
are controlled by the activities x and t of the horizontal and vertical edges of the 
lattice, respectively; the activity of the slanted edges in the circles, Fig. la, is set 
to unity. The statistical weight of each allowed set of trajectories is determined by 
the total number of horizontal and vertical unit steps in it. Thus, the generating 
(or partition) function of the model on a finite square lattice with L columns and 
M rows is given by: 

{C} 

Here the summation runs over the set of all allowed dimer configurations on the 
hexagonal lattice. For periodic boundary conditions, the solution of the Kasteleyn 
model in the thermodynamic limit is 



1 1 

/oo= lip \nZL^M{x,t) = — [ [dadp\n\l-te"'-xe'''\ (2) 
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The phase diagram of the model has four regions in the plane < t < oo, 
< X < oo, [T^\, denoted here by A,B,C and D. The phase in region A, 

{x,t) e A = {t> l+x} (3) 
describes stopped, close-packed cars. There is a trivial empty phase in region B, 

{x,t) e B = {x <l , t <\l- x\} (4) 

and the phase in C 

{x,t) G C = {a; > 1 , t < |1 (5) 

corresponds to the zero-density limit of cars moving at infinite speed. Normal tra- 
jectories of cars at all intermediate densities exist in the region {x,t) G D, 

{x,t) e D = {\l-x\ <t <\l+x\} (6) 

which motivates our interest only in that nontrivial case. 

The integration over a in equation (|]) can be easily performed by using the 
Jensen formula [Q. In region D the result is 

1 ^ 

/oo = A ln(t) + — [ dp ln(l - 2x cos(/3) + x^) (7) 
In J 



where 



, 1 fl+x^-t^' 

A = — arccos 

n \ 2x 
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Let us turn now to the interpretation of the observables in our model. For a 
finite lattice under periodic boundary conditions, the total numbers of steps in time 
and space can be expressed as 

Nt{C) = MNaiC) (9) 

and 

N,{C) = LNa{C)w{C) (10) 

where Na is the number of automobiles in configuration C, and w{C) is the number of 
windings of each trajectory in the periodic spatial direction. Therefore, the average 
density of cars can be defined as 

p(^, t) = YM ^^'^^^^ = ^WM""' = ^^""^ ^^^^ 

The rigorous definition of the average velocity v of cars is given by the average 
value of the ratio of space steps to time steps. However, as a first approximation we 
can assume 

. = (12) 

The numerator of the above ratio is obtained by differentiation of expression (|^) 
with respect to the activity of the space steps x\ 

^ (iV,(C)) = X j-/oo(x, t) = l{l- A(a:, t)) + ''^""^l K{x, t)) (13) 



LM ' ' " dx' ' ' ' 2' ' ' 2 
Here we have introduced the function 



r^/ X 2 I x+1 

K[x, tj = — arctan 



TC \ \X 



t2 - (1 -x)2 
-1|^ (l+x)2-t2 



(14) 



To obtain the fundamental traffic flow-density diagram, we start with the deflni- 
ton of the flow 

q = vp (15) 
which, in combination with (p!Tl)-(|T^), gives 

(iV.(C)) 

g(^,t) = -^ (16) 

The parameter t plays an auxiliary role in our model and we will express it in terms 
of the density p and activity x by using the relation A{x,t) = p, see (H) and ([TI|). 
Thus, for the traffic fiow q{x,p) we obtain the following explicit expression 



, , 1 — P siqnix 
q{x, p) = + - 



2n 



,(x- 1)2(1 + cos(7rp)) 

TT — arccos( — ; — 

x'^ — 2x cos(7rp) + 1 



(17) 
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Note that the variable x controls the average velocity Vm of a single car. Indeed, 
by considering configurations with only one trajectory, one can readily see that 
x/t. 

The fiow-density diagram at different fixed values of x is shown in figure 2. A 
remarkable feature of that diagram is the existence of two qualitatively different 
traffic-flow regimes. For x < 1 the flow reaches its maximum value 



(^(x, Pmax) 2 Pmax 



at the density 



Pmax = - arccos(x) (19) 

TT 

In the other regime, when x > 1, the traffic fiow decreases monotonically with the 
increase of p, from q{x, 0) = 1 to q{x, 1) = 0. 

The phenomenon of self-organized criticality in the transport fiow is clearly man- 
ifested in the asymptotic behaviour of the pair correlation functions both in space 
and time. Obviously, the correlations between time and space steps of the trajec- 
tories in our model correspond to dimer-dimer correlations in the original model of 
close-packed dimers on the hexagonal lattice. The latter can be obtained and studied 



by using the technique of Fisher and Stephenson |]12| . In the transport-fiow problem 
one is interested mainly in the correlations between cars passing a fixed position in 
space at two moments of time, say, t = and t = T, as well as in the correlations 
between cars at the same moment of time, but at a distance R apart. The former 
case corresponds to the correlation function for two vertical dimers in the same row 
of the decorated square lattice, and the latter case corresponds to two horizontal 
dimers in the same column, see Fig. 1. Due to the symmetry of the problem, we 
give here the expression for the temporal correlations only. By using the method of 



1^] , for the correlation function KtiT) one obtains in the thermodynamic limit 

27r 



^ ^ 



(20) 



The simple analytical structure of the integrand in the complex plane permits the 
exact evaluation of the above integral 

sin^fvrTA') , , 

Here A' = A'{x,t) = A{t,x). 

Similarily, for the spatial correlation function we obtain the asymptotic behaviour 
Ks{R) oc at large R. The power-law decay of both temporal and spatial corre- 
lation functions refiects the phenomenon of self-organized criticality in our model. 
However, the following questions arise: (1) are the obtained critical exponents uni- 
versal for the class of statistical models of traffic; (2) will they be the same for the 
cellular automaton models. 
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3 Comparison with the cellular automaton 
approach 

The interpretation of the above resuhs as a cellular automaton with local rules is not 
straightforward. The summation over all allowed dimer configurations corresponds 
to a global analysis of the statistics and there are no simple probabilistic local rules 
yielding the same distribution as the partition function (1). But it is possible to 
define a model with very similar properties as is shown in the following. 

Consider a probabilistic process according to the following rules. One time step 
consists an update of the positions of all cars in a parallel way. If a car has n empty 
sites in front of it (next car at n + 1) there are at most n partial movements of one 
site possible for this car. The decision how far the car finally moves depends on a 
series of random numbers yielding with probability p a '1' and with q = 1 — p a. '0'. 
The car moves as long as I's are chosen. When a occurs the movement stops until 
in the next time step this process starts anew. So if the first chosen random number 
is a the car does not move at all. If the series has I succeeding I's the car ends up 
at the l-th site (possibly directly behind the car ahead if it does not move itself in 
the same time step). 

A important feature of these rules is that no memory for the velocity is needed. 
On the other hand this can lead to unphysical features through enormous velocity 
fluctuations from one time step to the next. Nevertheless one can define a mean 
maximum velocity of a free driving car according to the rules just mentioned which 
is simply given hj v = p/q. This is the average distance a single car moves on a free 
road without other cars. 

In the following we show that the mean-field ansatz already gives the exact 
stationary state of the cellular automaton for any system size L and number of cars 
A^. The exactness of the mean-field solution has in this case the consequence that 
all possible states of the system are equally probable. 

We characterize a state of the system at time t only by the number di {i = 
1 , . . . , N) of free sites between car i and car i + 1 since no other quantities are 
needed. For the ring geometry this means that (In is the distance between car 
and car 1. It is then rather simple to write down the master-equation for the time 
evolution of the system 



P,+l({d,}) = E E ■ ■ ■ E n [P^^ (^'^.^.+1 + (1 - Sd.,s.^.)Q)] + 5.4-1}). 

Si=0s2=0 S]v=Oi=0 

(22) 

The variable Sj denotes the distance car i drives in the timestep t —>■ t + 1. The 
Kronecker-(5 takes into account the fact that if car i ends up directly behind the 
next car i + 1 (i.e. Sj+i = di) it does not have the 'deceleration' factor q. This is the 
main difference to the dimer-model where one has in any case this factor q (i.e. t in 
the other notation). 

The static solution of (^) is simply given by P{{di}) = i.e. asymptoti- 
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cally all probabilities are equal. This can be seen easily from 



1 = E E ■ ■ ■ E n [p^^ {^d.,s.,, + (1 - 5d.s.,Jg)] (23) 

si=0s2=0 S]v=Oj=0 

which follows from the identity 

l=j:p''{S,^..,s. + {l-S,,_,,sJq) (24) 

Si 

after interchanging product and sums. The probability /„ for a car to drive exactly 
n sites in the thermodynamic limit ^ oo is given by 

/„ = (l-p)V[(l-p)g + p] (25) 

since a car can stop by chance or since another car blocks it. For the flow in the CA 
one then obtains 

f \ ^ f P(l - P)P Pi^-P)v 
QcAiP, P)=P2^nfn = = _ 26 



To compare the flow (^) of the cellular automaton with the flow (|18]), we have 
to rescale (|T8|) with a factor 2 {q{x,p) 2q{x,p)). To illustrate this, we consider 
the limiting case of low car density and x —>■ 1. The dotted line in Fig. 2 tends to 
0.5 when p ^ 0. This reflects the fact that each horizontal line of the lattice can 
be occupied by the trajectory of the single car or be empty with equal probability. 
In the cellular automata approach,the motion with infinite velocity is forbidden and 
the case x — 1 corresponds to the maximal occupation of all horizontal bonds. This 
leads to the shift of the dotted line in Fig. 2 to the diagonal of the square. 

After the rescaling the maximum of both flows lie on the curve 1 — 2p. Therefore 
we can find the relation between the parameters p and x in the two models by 
identifying curves which have the same maximum. This yields the identification 

1 — - arccos x , , 

V=T. (27) 



(1-- 



arccos x 



Fig. 4 shows qcA{p, p) and q{x, p) for p = 3/4 and x = 1/2. Although the slopes for 
small and high densities are different the overall agreement of the two curves is very 
good. 

In conclusion we considered the Kasteleyn model as a model of traffic flow. This 
model allows to establish a bridge between the theory of exactly solvable models 
and the investigation of traffic flow by cellular automata. By choosing the activities 
in an appropriate way it is possible to interpret the dimer configurations as space- 
time trajectories of vehicles. The Kasteleyn model can be related to a modified 
Nagel-Schreckenberg model in which the cars do not have a velocity memory. The 
fundamental diagram of this modified CA can be obtained exactly as well and shows 
good agreement with the fundamental diagram obtained for the Kasteleyn model. 
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Figure Captions 



Figure 1: Representation of the hexagonal lattice as a decorated quadratic lattice, 
fig. la, and its mapping onto the square lattice, fig. lb. The encircled pairs of sites 
of the decorated lattice map on a single site of the square lattice, the vertical (hori- 
zontal) dimers on the former lattice map onto space (time) steps of the trajectories 
on the latter. 

Figure 2: The flow q{x,p) for different values of x. Curves with x < 1 lie below the 
dotted line (x = 1), curves with x > 1 above. The broken line indicates the location 
of the maximum flow. 

Figure 3: The flow qcAip, p) for different values of p. 

Figure 4: Comparison of ^^^(p, p) (dotted line) and 2q{x, p) forp = 3/4 and x = 1/2. 
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(a) (b) 

Figure 1: Representation of the hexagonal lattice as a decorated quadratic lattice, 
fig. la, and its mapping onto the square lattice, fig. lb. The encircled pairs of sites 
of the decorated lattice map on a single site of the square lattice, the vertical (hori- 
zontal) dimers on the former lattice map onto space (time) steps of the trajectories 
on the latter. 
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Figure 2: The flow q{x,p) for different values of x. Curves with a; < 1 he below the 
dotted line {x = 1), curves with x > 1 above. The broken line indicates the location 
of the maximum flow. 
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Figure 3: The flow qcA{p, p) for different values of p. 



l.Op— ^ 




p 



12 



Figure 4: Comparison oiqcA{p, p) (dotted line) and 2g(x, p) forp = 3/4 and x — 1/2. 




13 



